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In this paper we investigate deterministic diffusion in systems which are spatially extended in 
certain directions but are restricted in size and open in other directions, consequently particles 
can escape. We introduce besides the diffusion coefficient D on the chaotic repeller a coefficient 
D which measures the broadening of the distribution of trajectories during the transient chaotic 
motion. Both coefficients are explicitly computed for one-dimensional models, and they are found 
to be different in most cases. We show furthermore that a jump develops in both of the coefficients 
for most of the initial distributions when we approach the critical borderline where the escape rate 
equals the Liapunov exponent of a periodic orbit. 


I. INTRODUCTION 

Diffusion is one of the most common and most im¬ 
portant phenomena of nature. In physics it is impor¬ 
tant for nondeterministic systems, governed by noise, for 



transient diffusion in deterministic systems. By this we 
mean diffusion in a system spatially extended in certain 
directions which is however not closed in other directions 
and particles can escape. The diffusion is usually defined 
in the limit of infinitely long times. However, the intrin¬ 
sic feature of transient chaos (see for a review [p^ , pd]| ) is 
the limited time a system is close to a chaotic repeller. 
Only during this time the system shows the character¬ 
istics we are interested in. One can have in mind, for 
instance, a channel in a two-dimensional mesoscopic sys¬ 
tem containing Lorentz-type scatterers which is infinite 
in one direction but has a finite width in the transverse 
direction. Sooner or later all particles are scattered away 
from the channel and get lost if we assume free bound¬ 
aries. During their motion in the channel the particles 
are scattered and can exhibit diffusion in the extended 
direction. So the question arises in which way this diffu¬ 
sion should be characterized. 

The basic problem is how to define the diffusion coef¬ 
ficient. One possibility is to compute the diffusion coef¬ 
ficient D for particles on the repeller. That requires the 
computation of the diffusion coefficient for trajectories of 
those particles staying for ever in the channel. Follow¬ 
ing this line we can go to the limit of infinite time and 
thus can proceed as in the permanent chaos case. This 
formulation, however, can not provide a full description, 
since the probability is zero that a trajectory is on the 
repeller. A quite different view is to compute the actual 
properties in the neighborhood of the repeller, which is 
the main goal of the present paper. 



FIG. 1. An example for map F with translational symme¬ 
try and escape. The inset shows the corresponding reduced 
map /. 

To avoid complications with continuous time and 
higher dimensions we discuss here one-dimensional dis¬ 
crete maps defined on the real line. The extension of 
our formulation to more general cases is straightforward. 
The map Xt^i = F{Xt) has translational symmetry 
F"(a; -I- 1) = F(a;) -I- 1 and is not defined for a periodic 
set of intervals, called in the following windows, from 
which the trajectory escapes (Fig. 1). It is assumed that 
the map is chaotic in the sense that there exists no stable 
periodic orbit. Such maps can be described by reduced 
maps defined 'm.U = [0,1] Let us denote by X 

the box of the interval [N, N A phase point (called 
particle in the following) being after t iterations in box N 
with position Xt can either remain in that box or move 
to one of the nearby boxes or get lost. Which of these 
possibilities occurs depends only on fractional part {X} 
of X, since using the new variables N = [X],a; = {X} 
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we get a dynamical system of the two variables 

xt+i = f{xt), Nt+i = iVt + A(a;t), (1) 

which is equivalent to the map F. Here f{x) = {F(a::)} is 
the reduced map (cf. Fig. 1), which has a window where 
trajectories escape, and the integer part A(a;) = [F"(a;)] 
determines the steps among the boxes. We call a system 
with escape open, otherwise closed. The reduced maps 
with escape are generalizations of those introduced in 
and have been studied also in a different context 

Let us follow the reduced trajectories of particles in 
the interval U. At the beginning the distribution of the 
initial points Xq is according to an arbitrary density Pq 
in U. The change of P* after each step of iteration is 
induced by / alone and given by the Frobenius-Perron 
operator Pt+i(x) = CPt{x) = Y.z-.f(z)=xPt{z)/\f'{z)\. 
Repeated application of C leads to an average decrease 
of the density and to a change of its structure at the 
same time. Finally a quasistationary state is reached: 
the density normalized in U after each step has converged 
to P, the density of the conditionally invariant measure 
[p^ , whereas the density itself is decreasing after each 
iteration by a certain rate, related to the escape rate k: 

CP = XcP, |lnAc|=/c. (2) 



II. THE DIFFUSION COEFFICIENTS OF 
TRANSIENT CHAOS 

To get the diffusion coefficient for open systems we 
have to know the shift St of particles among the boxes 
at a given time 


t-i 

St=Nt-No = J2Mxk). (3) 

k=0 

Let us compute the average value of St first 
by averaging over an initial distribution Po(a:o) of 
xq values with P(xq) in U. Thus we obtain 
(5t)o = T!iZoIu,^ifix))Poix)dx/J^^Poix)dx with 

Ut = f~*(U) C\U. For t ^ Qo Ut becomes a repelling 
invariant set of the reduced map which will be called re- 
peller in the following. To get the diffusion coefficient 
the procedure is to take the limit t —> cx), which re¬ 
stricts the trajectories on the repeller. Apart from that 
the derivation of the corresponding expressions becomes 
straightforward along the lines of permanent diffusion 
[^,||. We use the Frobenius-Perron operator C, to write 
('5t)o = E‘=d/w[^*”'^'^'^o](a:)da://^[£*Po](a:)da;. But 
limj^oo acts as a projector on P, the density of the 

conditionally invariant measure. Thus the result for the 
averaged drift is 


This emphasizes the importance of the quasistationary 
state. Computing the properties of that state corre¬ 
sponds to taking the limit of large times T to reach it. 
We can define a coefficient D in the quasistationary state, 
which is analogous to a usual diffusion coefficient in the 
sense that it measures the mean quadratic deviation of 
the trajectories from their average shift. We call D the 
transient diffusion coefficient. In general D differs from 
D as will be shown explicitly for ID models. 

We study the properties of both diffusion coefficients 
as a function of a control parameter while keeping the 
escape rate constant. Within the family of such maps 
there is a limit when the Liapunov exponent of one of 
the unstable periodic orbits has reached the value of the 
escape rate and still the conditionally invariant measure 
has a density [^. This situation is called the critical 
borderline of transient chaos |H We show that here 


both D and D can have two different values depending on 
the initial distribution of trajectories. This will be shown 
to be traced back to the existence of two conditionally 
invariant densities with different escape rates. 

In Section 2 we derive the diffusion coefficient on a re¬ 
peller, D, and the transient diffusion coefficient D. Then 
we choose maps F of open systems whose reduced maps 
are in the family of fully developed chaotic maps . 

We present results that clearly demonstrate the differ¬ 
ence between D and D. In Section 3 we show that on 
approaching the critical borderline of transient chaos a 
jump develops in both D and D for most of the initial 
distributions. The conclusion ends the paper. 


,lim j(S't)o = / Mx)Kdx) = {A) (4) 

I—>00 I 

where fj,{dx) = 1 and fj, is the natural invariant mea¬ 
sure on the repeller. (S'Dq can be computed in an anal¬ 
ogous manner. The diffusion coefficient is then dehned 

by 


D = 


lim 

t—^OC) 


{{St-{St),?), 

2t 


( 5 ) 


which can be written in terms of the correlation function 
Caa on the repeller as 
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D = ^ CAA{\l-m\), ( 6 ) 

t^OO It ^^ 

/,m=0 

where C'aa(0 = {{A{f\x)) - (A)^)A(a;)^ . The ex¬ 
pressions are conmletely analogous to the formulae of 
closed systems [^,||. In spite of the analogies (A)^ and 
D cannot be accepted as the only quantities characteriz¬ 
ing open systems since a trajectory moves on the repeller 
with probability zero. 

Because of this a second procedure is suggested. Let 
us investigate particles that have been injected into the 
system and have survived (i.e. have not escaped) after 
a time T long enough to have the density Pt{x) of the 
reduced distribution of these particles close to the density 
P{x) of the conditionally invariant measure. We consider 
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the shift of the particles from time T until their escape 

Se{xT) = Nt+t{xt) - Nt = ^(^r+fe), where 

t{x) is the number of steps taken by the trajectory of 
cc, i.e. /’’(x) & U = [0,1] and ^ U. The 

average of this shift {Se)j' = jj^Se(x)PT{x) dx over the 
trajectories considered and the average lifetime (t)^ = 
Ik t(x)Pt(x) dx can be evaluated in the limit T —>■ oo 
using properties of the Frobenius-Perron operator. The 
result is (r) = Ac/(1 — Ac) and an average drift 

^ = A{x)P{x)dx. (7) 

v) Jui 

For sake of simplicity the oo index at the sign of average 
is omitted. 

We can characterize the transient spreading of the par¬ 
ticles by the transient diffusion coefficient defined as 


b = 


lim 

T—>oo 


2{t)t 


2(r)’ 


( 8 ) 


where Se^xr) = DILT iA{xT+k) - (A)), the cor¬ 
rected shift is introduced. Through steps similar to 
those which led to Eq. (|^) D can be related to the 
correlation function of transient chaos [|^ CAAin) = 
A-(n-i-i) j^^^^(A(jTt(^)) _ ^A))A(x)P(x) dx by the ex¬ 
pression 

CX> 

D= ^i”'cAA(|n|). (9) 


Physically this means the following. After injection drift 
and spreading of the injected particles depend on the ini¬ 
tial distribution, but its influence decreases exponentially 
fast. After a long time T it is no longer important and 
then further drift and spreading of the particles is dic¬ 
tated by the properties of the open system alone, charac¬ 
terized by (Se) and P, respectively. This general picture 
has to be refined, however, if the system has more than 
one smooth conditionally invariant measure as will be 
seen later. 

In the numerical calculations of open systems we take 
the reduced map / from the family of maps defined on the 
interval [0,1] and given by the inverse branches /;~^(x) = 
(x -I- v{x))/2R (lower branch), /^^(x) = 1 — /;~^(x) (up¬ 
per branch), where R > 1, x(0) = 0, v{x) is symmetric 
and |x'(x)| < 1. The density of a conditionally invariant 
measure can be expressed as 

P{x) = 1-1- v'{x). (10) 

The width of the window and the escape rate are given 
by 1 — R~^ and k = Ini?, respectively. Note that in the 
limit V = 0 the open tent map emerges whereas in the 
limit v'(0) = 1 we reach the state what we call critical 
borderline of transient chaos, since the Liapunov expo¬ 
nent of the fixed point at the origin equals the escape 
rate 00 • 


In our numerical computations we take ii = [0,1] and 
we use the piecewise parabolic map 00 which corre¬ 
sponds to v(x) = d ■ x(l — x), —1 < d < 1. We require, 
that a particle hitting the left (right) branch of / hops to 
the left (right) adjacent box N in the map F. If it hits 
neither branch then it escapes. This yields A(x) = —1 if 
X < l/(2i?) and 1 if x > 1 — l/(2i?). To be distinctively 
in the repeller regime we choose R = 1.5. However, the 
computation of D is not trivial in that regime. Reliable 
results were obtained only after we had applied back¬ 
ward iteration up to a depth of 20 steps, left out the 
corresponding gaps of the repeller and had started with 
10® trajectories (10® in case of d = 0.95). While com¬ 
puting b we did not take into account the first 20 steps 
(T = 20) to ensure the quasistationarity of the system. 
In Fig. 2 the diffusion coefficients D and b are shown 
for d values ranging from 0 to 1. Whereas they agree 
for the tent map a difference between D and b develops 
when the parameter d increases. The two extra points 
for d = 1 in Fig. 2 will be explained in the next section. 



FIG. 2. Numerical results for D and D in function of d. 
The superscript refers to the choice of the interval on which 
the initial points were distributed. 


III. DIFFUSION IN CRITICAL TRANSIENT 
CHAOS 

A peculiar situation emerges when we reach the crit¬ 
ical borderline of transient chaos (d = 1 in the exam¬ 
ple treated). Eq. (0 is still valid and gives a den¬ 
sity P^{x) (equal to 2(1 — x) in case of our example) 
for a conditional invariant measure with an escape rate 
= Ini?. The reason for introducing the superscript 
A is explained below. The corresponding eigenfunction 
P+^ of the adjoint Frobenius-Perron operator consists 
of two d-function peaks = (d(x) -I- d(l — x))/2 [^. 
This makes possible the existence of an other condition¬ 
ally invariant measure with density. To find it one should 
iterate a function (po that is positive everywhere in (0,1) 
and <(> 0 ( 0 ) = <(> 0 ( 1 ) = 0 by the Frobenius-Perron operator 
with the prescription (pt+i = atC(f)t, where at is chosen 
in such a way that the norm of 4>t+i is kept to unity. The 
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iterated functions remain orthogonal to because in 
the critical state the derivative of the map at the maxi¬ 
mum points is infinity. Furthermore it is easy to see that 
the positivity of is also maintained. In the limit t —> oo 
the iteration leads to a nonnegative eigenfunction which 
can be taken as the density P^{x) of a conditionally in¬ 
variant measure. The obtained eigenvalue defines a 
new escape rate hb = | In A^j. 

This means that we have two different conditionally 
invariant measures with different escape rates in the crit¬ 
ical state. Which of them is obtained in the long time 
limit when iterating the Frobenius-Perron operator de¬ 
pends on the initial density Po{x). It is clear from the 
above considerations that starting with a smooth Po{x) 
such that Po{0) = Po{l) = 0 we reach P^{x), otherwise 
P^{x). When simulating the initial density by a large 
number of initial points the obtained distribution Pt{x) 
of iterates at a large time T is governed by the same rule. 
Note that to balance the escape we have to increase the 
number of initial points when increasing T. 

The above rule for the asymptotic distribution is re¬ 
flected in the values of the diffusion coefficients D and D. 
In the numerical calculations we distributed a large num¬ 
ber of particles uniformly on intervals, [0,1] and [0.1, 0.2], 
to simulate initial densities in the basin of attraction of 
P"^ and P^, respectively. It is clear from Fig. 2 (at d = 1) 
that the result depends strongly on the initial distribu¬ 
tion (as well as values of the escape rate and the average 
drifts (A)^, (A), not displayed here). Note, the value 
of is zero due to the fact that the natural measure 
induced by P^ on the repeller (but not the repeller it¬ 
self) degenerates, it becomes a Dirac delta function at 
the origin 011 ■ 

IV. DISCUSSION 

In this paper we have investigated statistical proper¬ 
ties in long time limit of trajectories in the transiently 
chaotic regime by the help of two diffusion coefficients D 
and D. For D this time limit refers to the whole length 
of the trajectories (a procedure leading automatically to 
experience the properties of the repeller itself). In con¬ 
trast the long time limit of D is only needed to measure 
it in the quasistationary state. Though the definition of a 
transient diffusion constant as introduced in the present 
paper seems to be the most natural one, it is not the 
only possibility. In future work we plan to examine also 
alternative definitions. 

We have shown furthermore that a map at the critical 
borderline of transient chaos has two smooth condition¬ 
ally invariant measures. We find, however, that in nearly 
all practical situations (i.e. when the initial phase points 
are concentrated in the inner part of the interval) the 
relevant measure is not the usual one, belonging to the 
smaller escape rate, but an other eigenfunction of the 
Frobenius-Perron operator C. 


In this paper we have treated a ID model which can 
be conceived as describing the dynamics along the ex¬ 
panding direction of a 2D map. It is expected that a 
more complete treatment would not change the conclu¬ 
sions. We plan as a next step to study 2D Hamiltonian 
systems. 
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